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ABSTRACT 

We propose a magnetic mechanism for the collimated explosion of massive stars relevant for long-duration 
gamma-ray bursts (GRBs), X-ray Flashes (XRFs) and asymmetric core collapse supernovae. In particular, we 
apply Lynden-Bell's magnetic tower scenario to the interior of a massive rotating star after the core has col- 
lapsed to form a collapsar with a black hole accretion disk or a millisecond magnetar as the central engine. 
The key element of the model is that the pressure of the toroidal magnetic field, continuously generated by 
differential rotation of the central engine, drives a rapid expansion which becomes vertically collimated after 
lateral force balance with the surrounding gas pressure is reached. The collimation naturally occurs because 
hoop stress concentrates magnetic field toward the rotation axis and inhibits lateral expansion without affecting 
vertical expansion. This leads to the growth of a self-collimated magnetic structure which Lynden-Bell termed 
a magnetic tower. When embedded in a massive star, the supersonic expansion of the tower drives a strong 
bow shock behind which an over-pressured cocoon of shocked stellar material forms, as observed in hydrody- 
namical simulations of collapsar jets. The cocoon confines the tower by supplying collimating pressure support 
and provides stabilization against disruption due to magnetohydrodynamical instabilities. Because the tower 
consists of closed field lines starting and ending on the central engine, mixing of baryons from the cocoon 
into the tower is suppressed. Baryon loading due to magneto-centrifugal winds from the central engine may 
also be suppressed because of the expected field line geometry. The channel cleared by the growing tower is 
thus plausibly free of baryons and allows the escape of magnetic energy from the central engine through the 
star. While propagating down the stellar density gradient, the expansion of the tower accelerates and becomes 
relativistic. At some point during the expansion fast collisionless reconnection becomes possible. The resulting 
dissipation of magnetic energy into accelerated particles may be responisble for GRB prompt emission. 

Subject headings: accretion, accretion disks — magnetic fields — MHD — gamma rays: bursts — supernovae: 
general — stars: magnetic fields 



1. INTRODUCTION 

Long duration gamma-ray bursts (GRBs) are asymmetric 
explosions (Harrison et al. 1999; Stanek et al. 1999) as- 
sociated with the death of massive stars (Matheson et al. 
2003; Hjorth et al. 2003; Stanek et al. 2003). Spectro- 
scopic observations of the afterglow of GRB030329 revealed 
a bright underlying Type Ibc supernova, SN2003dh (Mathe- 
son et al. 2003; Stanek et al. 2003), which bore close re- 
semblance to SN1998bw, the first supernova associated with 
a GRB, the faint nearby (40 Mpc) GRB980425 (Galama et 
al. 1998). Two other convincing spectroscopic identifications 
of supernovae associated with XRF020903 (Soderberg 2005) 
and GRB031203 (Malesani et al. 2004) have been reported 
in addition to about ten photometric detections (Zeh, Klose 
& Hartmann 2004; Soderberg et al. 2006a). It is now firmly 
established observationally that most, perhaps all, of the long 
duration GRBs are associated with Type Ibc supernovae. 3 Ev- 
idence for strong asymmetry in GRB ejecta is derived from 
the afterglow light curves which exhibit, after several days, an 
achromatic transition in the power law decay of brightness as 
a function of time, indicating that the explosion was beamed 
into jets with opening angle of ~ 10 degrees. Taken together, 
these observations show that some supernova explosions in- 
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volve strongly asymmetric relativistic outflow. Massive stars 
are thus observed to be capable, in some cases, of producing 
strongly collimated energetic outflows when they die. 

Such asymmetry is naturally expected in models where the 
progenitor star is rapidly rotating when it collapses and sub- 
sequently explodes. Stellar rotation breaks spherical sym- 
metry and provides the rotation axis as a preferred direc- 
tion along which jet-like outflow can develop. Observations 
of massive stars do indicate rapid surface rotation velocities 
(Fukuda 1982), though the distribution of angular momentum 
in the cores of evolved stars when they undergo core collapse 
is uncertain (Spruit 2002; Heger, Woosley & Spruit 2005; 
Maeder & Meynet 2005). However, low metallicity is ex- 
pected to be beneficial for maintaining angular momentum in 
the core until collapse (MacFadyen & Woosley 1999; Hirschi, 
Meynet, & Maeder 2005). Current stellar evolution models 
do, in some cases, predict rapid rotation in collapsing cores 
of massive single stars (Heger, Woosley & Spruit 2005), per- 
haps aided by global mixing (Woosley & Heger 2005; Yoon 
& Langer 2005). It is also expected that membership in a 
tight binary system will endow some massive stars with rapid 
rotation in their cores when they die (Fryer & Heger 2005; 
Petrovic et al. 2005) though the angular momentum imparted 
to the core depends on detailed transport processes (Petro- 
vic et al. 2005). 

A theoretical model for long duration GRBs based on rapid 
stellar rotation is known as the collapsar (Woosley 1993; 
Paczynski 1998; MacFadyen & Woosley 1999; MacFadyen, 
Woosley & Heger 2001) in which the core of a massive rotat- 



2 



ing star collapses to form a black hole. Numerical simulations 
have shown that, in rapidly rotating stars, an accretion disk 
rapidly forms around the young black hole (MacFadyen & 
Woosley 1999). Because of angular velocity gradients in the 
disk, the magneto-rotational instability (MRI) is expected to 
develop, providing angular momentum transport and dissipa- 
tion of gravitational energy. At the temperatures and densities 
(r~4x 10 10 K, p ~ 10 g cm -3 ) present in collapsar disks, 
neutrino emission can cool the accreting gas with a range of 
efficiency. The collapsing outer stellar core is therefore able 
to accrete at rates of ~ 0.1Mqs~ 1 for times > 10 s provid- 
ing sufficient energy for sufficiently long times to power long 
GRBs. 

Magnetohydrodynamical (MHD) simulations have demon- 
strated the development of MRI in collapsar disks (Proga et al. 
2003; McKinney 2005). Magnetic fields approaching 
equipartition values in excess of 10 15 G are then expected 
to develop as the disk forms. In addition, field strengths 
approaching 10 15 G may develop during collapse as seed 
poloidal field winds up into toroidal field. Similar pro- 
cesses, resulting in similar field strengths, are also believed 
to be taking place in core-collapse supernovae explosions 
(Akiyama et al. 2003; Ardeljan, Bisnovatyi-Kogan, & Moi- 
seenko 2005). 

In addition, it is possible that a turbulent dynamo is not re- 
quired in the collapsar case. This is because a plausible source 
of large-scale poloidal field threading the disk may be the field 
advected and compressed from the pre-collapse star. Mag- 
netized white dwarfs are observed to have surface magnetic 
fields of 10 9 G. If the cores of some massive cores posses 
fields of similar magnitude, they will grow to 10 14 G if a core 
of radius 10 9 cm and B ~ 10 9 G collapses to 3 • 10 6 cm. As 
was recently established by Braithwaite & Spruit (2005), sta- 
ble large-scale poloidal field structures can be maintained for 
a stellar lifetime by a spheromak-like field configuration in- 
volving twisted toroidal field inside the star. 

Because gas along the rotation axis of the star does not 
experience a centrifugal barrier, it falls into the black hole 
unimpeded, thus lowering the density above the hole to p < 
10 6 gm cm" 3 . Similarly to its role in Type II supernovae ex- 
plosions (Goodman et al. 1987), neutrino annihilation may be 
important in the low density funnel above the black hole, heat- 
ing the stellar gas and causing it to expand along the rotation 
axis of the star. Under favorable conditions this hot gas can 
form a collimated relativistic outflow which is capable of es- 
caping the stellar surface and forming a fireball. MacFadyen 
& Woosley (1999) explored this possibility but noted that 
magnetic processes are an alternate mechanism capable of ex- 
tracting energy from the system. In fact, the suggestion that 
strongly magnetized jets play an active role in GRB explo- 
sions (perhaps via a version of the Blandford-Znajek mech- 
anism) has been invoked by several authors (e.g., Thomp- 
son 1994; Meszaros & Rees 1997; Lee et al. 2000; Vla- 
hakis & Konigl 2001 ; van Putten & Ostriker 2001 ; Drenkhahn 
& Spruit 2002; Lyutikov & Blandford 2003; van Putten & 
Levinson 2003; Proga et al. 2003; Lyutikov 2004, 2006; 
Lei et al. 2005; Proga & Zhang 2006). In addition, a lot 
of relevant work has been done in the supernova context, 
where several magnetic explosion mechanisms have been ex- 
plored (e.g., Leblanc & Wilson 1970; Meier et al. 1976; 
Wheeler et al. 2000, 2002; Akiyama et al. 2003; Ardel- 
jan et al. 2005; Blackman et al. 2006). 

In this paper, motivated by the above considerations, we 
consider a possible magnetic mechanism for driving a low- 



baryon-load, Poynting-flux dominated jet through a massive 
collapsing star. In order to understand the structure of mag- 
netic jets in massive stars we investigate a simple analytic 
model that may help guide interpretation of detailed MHD 
simulations. 

Our paper is organized as follows. In §|2]we describe the 
key physical ideas that underlie our model. The main idea 
is that we apply the magnetic tower mechanism, developed 
previously by Lynden-Bell (1996) for active galactic nucleus 
(AGN) jets, to the collapsar environment. We first give a qual- 
itative description of Lynden-Bell's original model (in an ac- 
cretion disk g eom etry) and derive some basic scalings result- 
ing from it (§ 12. 11 1. In the rest of §0 we present our picture 
of a magnetic tower growing rapidly insid e a s tar and driv- 
ing a strong shock through it. Thus, in § 12.21 we describe 
our modifications to the original magnetic tower model; the 
most important of these is that the external pressure confin- 
ing the tower is no longer arbitrary but is determined self- 
consistently by the hot gas cocoon surroundin g the tower and 
by the strong shock driven into the star. In § 12.31 w e der ive 
simple scalings for the tower parameters and then, in § 12.41 we 
use these scalings to make some quantitative estimates. In §|3] 
we make the next step and present mathematical formalism of 
our model. Our analysis is characterized by representing the 
time evolution of the tower by a sequence of force-free mag- 
netostatic equilibria, obtained by solving the Grad-Shafranov 
equation with additional constraints and boundary conditions 
that, in general, change with time. We then present, as exam- 
ples, several particular analytical solutions (in §0}. After that, 
in § |5J we discuss the implications of our model for gamma- 
ray bursts and outline some important open issues that we feel 
need to be addressed in future research. In particular, in § 15.11 
we argue that, as the tower expands in a stratified star, the un- 
perturbed stellar density at the top of the tower drops and the 
expansion process accelerates. At a certain point, the tower 
growth will then have to make a transition to the relativistic 
regime, which is not covered by our theory. Nevertheless, we 
speculate that the tower will remai n col limated, with an open- 
ing angle of a few degrees. In § 15.21 we discuss the tower 
structure a nd th e distribution of magnetic energy inside the 
tower. In § !5.3l we discuss the issue of baryon contamination 
originating from the disk as a result of ne utrino ablation and 
magneto-centrifugal winds. Next, in § 15.41 we discuss various 
MHD instabilities that may be excited in the system and the 
effe ct th ey may have on the evolution of the magnetic tower. 
In § 15.51 we discuss the question of whether the tower can be 
disrupted through reconnection while it is still d eep inside the 
star and argue against this possibility. In § !5.6l we outline nu- 
merical simulations which may be performed to explore the 
model. We present our conclusions in §[6] 

Finally, note that for definiteness we develop our model in 
the context of the original collapsar scenario, where the cen- 
tral engine is an accretion disk around a black hole. The en- 
ergy source in this case is accretion energy. However, we be- 
lieve that a similar magnetic tower mechanism may also work 
when the central engine is a young millisecond magnetar born 
inside a collapsing star. The typical values of the magnetic 
field strength, the rotation rate, and the size of the system are 
similar in the two cases. The overall electro-magnetic lumi- 
nosities should therefore also be comparable (e.g., Usov 1992; 
Thompson 1994). The energy source for the explosion in the 
magnetar case is the rotational energy of the neutron star and 
the magnetic field acts as an agent that extracts this energy. 
The basic mechanism is essentially similar to that proposed by 
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Ostriker & Gunn ( 1 97 1 ) for powering supernovae light-curves 
by the spin-down magnetic luminosity of a rapidly-rotating 
pulsar (with the magnetic field scaled up by three orders of 
magnitude and the timescale scaled down by six). In fact, sev- 
eral models invoking rapidly-rotating magnetars as central en- 
gines for GRBs have been proposed (e.g., Usov 1992; Thomp- 
son 1994; Yi & Blackman 1998; Thompson et al. 2004; 
Lyutikov 2006). However, these models usually don't dis- 
cuss the geometry (e.g., collimation) of the outflow. In ad- 
dition, they typically consider an isolated magnetar (without 
a surrounding stellar envelope). The collimating influence of 
the dense stellar material is thereby neglected. We believe 
that a rapidly -rotating magnetar will make an even better cen- 
tral engine when considered within the collapsar framework. 
A simplified physical model representing this scenario is the 
pulsar-in-a-cavity problem. By analyzing this problem, we 
can show that, even though the magnetar itself rotates uni- 
formly, a strong differential rotation is effectively established 
on the field lines extending beyond the light cylinder. As a re- 
sult, toroidal magnetic flux is constantly injected into the cav- 
ity and eventually toroidal magnetic field becomes dominant 
over both the poloidal magnetic field and the electric field in 
the magneto sphere. Any subsequent expansion of the cavity 
is then going to be mostly vertical because of the collimating 
hoop stress, just as in Lynden-Bell's (1996) model. Even- 
tually, a magnetic tower forms. In our opinion, this pulsar- 
in-a-cavity problem, viewed as a paradigm for a millisecond- 
magnetar-driven GRB explosion mechanism, is very impor- 
tant by itself. Therefore we believe it deserves a separate 
study and we intend to pursue it in an upcoming paper. 

2. MAGNETIC TOWER: BASIC PHYSICAL PICTURE 

2.1. Lynden-Bell's Original Magnetic Tower Model 

As a specific mathematical vehicle to illustrate our ideas, 
we choose the magnetic tower introduced by Lynden-Bell 
(1996, 2003). A magnetic tower is an axisymmetric magnetic 
configuration that arises when a system of nested closed flux 
surfaces, anchored in a differentially-rotating disk, is twisted 
and, as a result, inflates, but when this inflation is controlled 
by a surrounding external pressure. The basic physical mech- 
anism of this process can be described as follows (see Fig.Q. 

Consider a thin conducting disk with some vertical mag- 
netic flux frozen into it. Let us assume that initially the mag- 
netic field is potential and has a fully-closed dipole-like topol- 
ogy (see Fig. [^), with the two footpoints of each field line 
located at different radii on the disk. Now let us assume that 
the disk rotates non-uniformly, which is certainly the case for 
a Keplerian disk. Then, each field line marked by the parame- 
ter \& is twisted at a rate Af2(\&) equal to the difference in an- 
gular velocities of the two footpoints of this line. Correspond- 
ingly, toroidal magnetic flux is generated from the poloidal 
flux. The pressure of this toroidal field pushes the flux sur- 
faces outward, against the poloidal field tension. It is assumed 
that during the initial stages of this process the gas pressure, 
as well as the gravitational and inertial forces, are negligi- 
bly small in the disk magnetosphere, so the magnetic field 
is force-free. Then the expansion is uncollimated, typically 
along the direction making a 60° angle with respect to the ro- 
tation axis (Fig.^i). However, as was shown by Lynden-Bell 
(1996), if there is some, no matter how small, external gas 
pressure that surrounds the expanding disk magnetosphere, 
then at some point the sideways expansion is stopped. This is 
because the magnetic field strength at the outermost portions 
of inflating field lines decreases rapidly during the inflation 
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FIG. 1 . — Development of a magnetic tower in Lynden-Bell's (1996) model. 



process; then, once B 1 /87r drops to the level of the external 
gas pressure, any further horizontal expansion ceases. How- 
ever, as again was shown by Lynden-Bell (1996), the story 
doesn't end here; unable to expand sideways, the twisted mag- 
netosphere expands in the vertical direction (Fig. QJ), even- 
tually forming a slender cylindrical column that Lynden-Bell 
termed a magnetic tower (see Fig. 01). If the external pressure 
outside of the tower is kept constant and uniform, then the top 
of the tower rises at a constant speed. Plasma inertia never 
plays any role in this process; the entire evolution is viewed 
as a sequence of magnetostatic equilibria, with the field being 
force-free inside the tower and with pressure balance between 
the external gas outside of the tower and the magnetic field 
inside. 

Note that the assumption that both ends of a field line con- 
nect to the disk itself is actually not essential. An alterna- 
tive configuration, with a very similar overall behavior, is 
that of a rotating conducting disk magnetically connected to a 
highly-conducting rotating central star (a proto-neutron star 
in the collapsar context) or even a rotating black hole. In 
the conducting star case, the differential star-disk rotation 
leads to the same field-line inflation and opening process in 
the force-free regime (e.g., Lovelace et al. 1995; Uzden- 
sky et al. 2002), followed by the tower stage when the exter- 
nal pressure becomes important. A very similar process takes 
place when a conducting disk is connected magnetically to a 
rapidly-rotating black hole: even though a black hole cannot 
be considered a good conductor, general-relativistic frame- 
dragging ensures that a closed-field configuration that extends 
far enough on the disk cannot stay in equilibrium and there- 
fore has to inflate and open up (Uzdensky 2005). This picture 
is consistent with that observed in recent general-relativistic 
MHD (GRMHD) numerical simulations of accreting black 
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holes (McKinney 2005; Hawley & Krolik 2006). 

In order to get accurate quantitative results, one needs to 
consider a specific solution of the governing MHD equations 
with some particular boundary conditions. However, to get 
a basic physical feeling of how the magnetic tower grows, it 
is instructive to derive some simple order-of-magnitude es- 
timates and scaling relationships. In Lynden-Bell's (1996) 
model, the main input parameters that determine the solu- 
tion are the total poloidal magnetic flux \&o (P er un it toroidal 
angle) in the tower, the characteristic differential rotation 
rate A£l, and the external pressure P ext . In the rest of this sec- 
tion we investigate how the main parameters that characterize 
the tower, namely, its radius, the typical magnetic field, and 
the growth velocity, scale with these three input parameters. 

First, we estimate the radius of the tower, Rq, and the char- 
acteristic poloidal magnetic field strength, fi po i. They are re- 
lated via 



B, 



pol 



fin 



(i) 



As we mentioned earlier, the radius adjusts so that the mag- 
netic pressure inside the tower equals the external gas pres- 
sure. From the force-free balance inside the tower we ex- 
pect the toroidal magnetic field to be roughly the same as 
the poloidal field: fi^ ~ fi po i ~ fin. The total magnetic field 
strength at the outer edge of the tower is thus also of the order 
of fin. Then, from the the condition of pressure balance across 
the tower's side wall we get 



fi ~ VSTTfiexf 

By combining this with equation Q, we get 



8vrfi e! , 



(2) 



(3) 



Next, let us estimate the rate of growth of the tower. To 
do this, we use the fact that the toroidal magnetic flux is con- 
tinuously generated out of the poloidal flux by the differential 
rotation. Every time a given flux tube with poloidal flux is 
twisted by one full turn, the amount of toroidal flux carried by 
the tube is increased by A "J. Therefore, after N = A£lf / (2tt) 
turns the total toroidal flux \ m the tower becomes 



X = 27r* N = * A£lf . 



(4) 



Assuming that this toroidal flux fills the tower uniformly 
and that the tower is a cylinder with radius fin and height Z top , 
we see that the characteristic toroidal field has to be of the 
order of 



finZtop 
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Aflt = fin — Ant . 
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However, as we stated earlier, the typical toroidal field in 
the tower is of the order fin; therefore, the height of the tower 
has to increase steadily as 

Z,op(f) ~ fi A£lf . (6) 

In other words, the tower grows at the speed comparable to 
the typical differential rotation velocity Ro An. If the external 
pressure does not change, the radius of the tower, determined 
by equation Q, stays constant during its growth; therefore, 
after many turns (Ant ^S> 1), the height Z top becomes much 
larger than the radius, i.e., the tower becomes slender. 

Since the first analytical solution proposed by Lynden-Bell 
(1996), the magnetic tower concept is becoming more and 



more accepted by the astrophysical community. For exam- 
ple, the formation and evolution of magnetic towers have been 
studied in numerical simulations (Li et al. 2001; Kato et al. 
2004) and even in real laboratory experiments (Hsu & Bellan 
2002; Lebedev et al. 2005). When comparing with the recent 
GRMHD simulations of accreting black holes (e.g., McKin- 
ney 2005; Hawley & Krolik 2006), one should note that there 
is no confining external gas pressure in these simulations, and 
so one may not expect outflow collimation by the magnetic 
tower mechanism. However, one may expect a poorly colli- 
mated tower-like magnetic structure on very large scales, with 
a narrow Blandford-Znajek (1977) jet consitituting the inner 
core of the tower. 

Energy Flow in Magnetic Towers 

An interesting question is how energy flows through a mag- 
netic tower. As one can easily see, Poynting flux flows up 
from the disk along the inner segment of each field line and 
down to the disk along the outer segments. In our view, this 
can be understood as follows. For each field line ^, the inner, 
faster-rotating footpoint (1) performs work on the magnetic 
field at a rate proportional to W\ ~ 7 po i[5'(l)]il(l) (per unit of 
poloidal flux). The corresponding decelerating torque per unit 
flux is proportional to t\ ~ /poit^Cl)]- In turn, the magnetic 
field exerts an accelerating torque per unit flux T2 ~ / po i[ v I'(2)] 
on the outer disk footpoint of the same field line. Correspond- 
ingly, it performs work at a rate W2 ~ 7 po i[ , I'(2)]£l(2). Because 
of force-free equilibrium in the tower, I po \[^>(l)] = I po i[^(2)] = 
Ipoii^), and so t\ = t%, i.e., all the angular momentum ex- 
tracted magnetically from point 1 is transferred to point 2. 
The two energy flows, on the other hand, are not equal. In- 
deed, since £1(1) > £1(2), the energy that is extracted from 
point 1 and that flows up along the inner segment of the field 
line is greater than the energy that flows down along the outer 
segment and is deposited in the disk at point 2. The differ- 
ence, proportional to I po i(^>), is the actual power that drives 
the expansion of the tower. A part of it goes into filling the 
growing volume of the tower with magnetic energy, and the 
rest goes into performing work against external gas pressure 
and driving the shock through the star. 

It is interesting to note that the total vertical Poynting flux in 
the two segments only involves the differential rotation An = 
n{\)- £1(2), but is independent of the absolute rotation itself. 
This is because we are dealing here with a force-free equi- 
librium, so that /poi is constant along the entire length of a 
field line; in particular, it has the same sign on the two seg- 
ments of the field line, and hence so does B^. The situation 
would be drastically different in the relativistic-rotation case, 
where both field-line segments extend beyond their respec- 
tive light cylinders. In that case, one would no longer have a 
force-free equilibrium along the entire stretch of a field line; 
in particular, equilibrium will most likely break down at the 
farthermost tip of the line where the two segments join. As a 
result, the signs of 7 po i (and hence of fi^) on the two segments 
would be opposite, which corresponds to both segments be- 
ing swept back. Consequently, the Poynting flux would be 
outward along both segments. A similar situation arises in the 
non-force-free MHD case; the two field-line segments would 
then be swept back by plasma inertia if they extend beyond the 
Alfven point. This would again result in an outward Poynting 
flux along both segments. In both of these cases, the total ver- 
tical Poynting flux depends on the absolute rotation rates £1(1) 
and £1(2) themselves, as opposed to just their difference. 



2.2. Magnetic Tower Driving a Shock through a Star 

As we mentioned briefly in the Introduction, we believe that 
there are several good reasons that make the magnetic tower 
mechanism an attractive model for the formation and propa- 
gation of a magnetically-dominated jet through a star within 
the collapsar model for GRBs (and perhaps core-collapse su- 
pernovae as well). To reiterate, a configuration where all the 
field lines close back onto the central engine (e.g., an accre- 
tion disk or a magnetar) is very natural for a field created by a 
turbulent dynamo (with zero net flux). In addition, diffusion 
of particles across the field is suppressed and thus baryon con- 
tamination of the tower from the stellar envelope is inhibited. 
This is in contrast with the models that invoke the Blandford- 
Znajek mechanism operating directly along open field lines 
inside the star. Even if a light, ultra-relativistic outflow is 
launched along such field lines, it will first have to push all the 
baryonic matter, already present on these field lines, ahead of 
itself; therefore, the resulting ejecta in these models will not 
be baryon-free. On the other hand, the closed-field geometry, 
characteristic of the magnetic tower model, at least eliminates 
this problem. The only way baryons can get into the tower 
is via a wind blowing from the disk or the neutron star. An 
assessment of baryon contamination due to a wind is an im- 
portant issue that we will discuss in § 15.31 

In order to apply the magnetic tower model to the collapsar 
scenario, we first have to make a few modifications to Lynden- 
Bell's original picture. Specifically, we extend his model by 
taking into account the high-pressure cocoon that surrounds 
and confines the tower. In our model, the magnetic tower 
grows very rapidly and thus acts as a piston driving a shock 
ahead of it. The shocked stellar gas above the tower has very 
high pressure; it therefore squirts sideways and forms back- 
flows that fill the cocoon around the tower as described for 
AGN jets by Begelman, Blandford & Rees (1984) and in the 
collapsar case by Ramirez-Ruiz, Celotti & Rees (2002) and 
Matzner (2003) (see Fig.|2jl. Therefore, the external pressure 
that confines the tower is no longer an arbitrary parameter, in 
contrast with Lynden-Bell's model. Instead, it is determined 
by the jump conditions across the shock surrounding the co- 
coon across the contact discontinuity between the cocoon and 
the tower. Since the speed at which the magnetic tower plows 
through the star is much higher than the sound speed in the un- 
perturbed stellar material, the external unperturbed pressure 
of the star is irrelevant; it should thus be excluded from our list 
of the three main input parameters used for our simple, order- 
of-magnitude estimates. Instead, the expansion of the cocoon 
is controlled by the ram pressure related to the gas inertia; 
therefore, we replace the external pressure by the unperturbed 
stellar density po in the list of basic dimensional parameters 
(along with Wo and A51) that determine the relevant physical 
scales in our problem. This change represents an important 
difference between our model and Lynden-Bell's. 

The actual situation is complicated further by the two- 
dimensional character of the problem. Indeed, since the sound 
speed inside the cocoon is very high, comparable to the speed 
at which the tower grows, gas pressure tends to be partially 
equalized throughout the cocoon. This expectation is sup- 
ported by the hydrodynamic numerical simulations 4 of the 
collapsar model by MacFadyen et al. (2001) and Zhang et al. 

4 Note that, even though the simulations by MacFadyen et al. (2001) and 
Zhang et al. (2003) do not include magnetic fields, they probably still give 
a reasonably good qualitative representation of the cocoon behavior even in 
our magnetic tower scenario, since the cocoon itself remains unmagnetized. 
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FIG. 2. — The main components of our model. The lines attached to the 
disk represent poloidal magnetic field lines (toroidal magnetic field is also 
present). A magnetic tower grows rapidly and drives a strong shock through 
the star. The shocked stellar gas behind the shock forms a hot cocoon whose 
high pressure confines the tower. 



(2003), which show that, typically, the variation of the cocoon 
pressure along its length is relatively weak, just a factor of 5 
or 10. This is very moderate compared to the correspond- 
ing variations of the unperturbed stellar density and pressure, 
which both vary by many orders of magnitude along the ver- 
tical extent of the cocoon. As a result, the gas pressure is 
very high everywhere in the cocoon so the cocoon also drives 
a sideways shock into the star. Thus, the boundary between 
the tower and the cocoon is a contact discontinuity, whereas 
the boundary between the cocoon and the rest of the star is a 
two-dimensional strong shock of some complicated shape. 

We see that, in principle, the problem calls for a considera- 
tion of the entire two-dimensional shock-cocoon-tower struc- 
ture. For simplicity, however, we shall represent the cocoon 
by only two components: the hot spot right above the tower 
(with some uniform pressure P top ) and the cylindrical shell 
surrounding the tower on the sides (with a different uniform 
pressure P s id e )- Thus, we shall parameterize the pressure non- 
uniformity in the cocoon by the ratio r\ = P s ide/Aop- In ac- 
cordance with this 2-component cocoon structure, we shall 
represent the strong shock between the cocoon and the star 
also by two components: a plane shock that propagates purely 
vertically above the tower and a slower cylindrical shock that 
propagates sideways. An important difference between these 
two shocks is that, whereas the upward shock is driven exclu- 
sively by the vertical growth of the magnetic tower, the side- 
ways shock is driven mostly by the pressure of the hot gas that 
flows into the side part of the cocoon from the top hot spot; 
the horizontal expansion of the tower may also play a role in 
driving the sideways shock, but probably to a lesser extent. 

Let us now discuss what determines ?/. We use sim- 
ple energetics considerations for a rough estimate. The en- 
ergy flux that is required to drive the sideways shock is of 
order P s ideVs.side, where V s , s id e is the speed of the sideways 
shock; it is, in turn, proportional to the sound speed, V s , s jd e ~ 
(Aide/ Po,side) (1,/2) - In our simple 2-component model, the area 
of the sideways shock is larger than the area of the vertical 
shock by a factor k = Z top 7? cocoon /7?Q ^> 1. Thus, in order for 
the power necessary to maintain the sideways shock not to 
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exceed the available thermal power generated by the vertical 
growth of the tower, one must satisfy 



Tj < K 



2/3 ( PO.side 
P(),top 



1/3 



(7) 



Thus, we see that, as the tower grows (and k increases) 
with time, 77 = P S ide/Ptop tends to decrease, but, on the other 
hand, this decrease is partially weakened by the density non- 
uniformity, i.e., by the large ratio po, s ide/po,top- We also see 
that the sideways shock propagates noticeably slower than the 
vertical shock. This is due to a combination of two factors: a 
smaller (by a factor of iff pressure driving the sideways shock 
and a higher background density of the stellar gas into which 
the sideways shock is propagating. As a result, the cocoon is 
moderately elongated in the vertical direction, with an aspect 
ratio of 



Zcocoon V^top — 1/2 / P0,side 



R 



V, 



s.side 



Pfltop 



P0,side 
P0,top 



1/3 



(8) 



These estimates, of course, are only as good as the under- 
lying two-component cocoon structure assumed here. In a 
realistic situation, where the tower is growing in a strongly 
stratified medium, the radii of both the tower and the cocoon 
are not going to be constant in the vertical direction. In that 
situation, one cannot represent the tower and the cocoon by 
straight cylinders; instead, one has to describe their shapes by 
some increasing functions Ro(z) and R CO coon(z)- 

2.3. Simple Estimates 

Let us now show how the basic properties of the growing 
magnetic tower are expressed in terms of the three parame- 
ters \I>, Afl, and po. For simplicity, in this and the next sub- 
sections, we ignore the non-uniformity of the gas pressure in 
the cocoon, i.e., we set 77 = 1. We will reinstate this non- 
uniformity and study the dependence of the tower growth pa- 
rameters on r) in our more rigorous analysis of §§|3j|4] 

For simplicity we shall use shock jump conditions corre- 
sponding to a one-dimensional problem. Since the gas pres- 
sure of the unperturbed stellar material upstream of the shock 
is neglected, the shock is strong. Assuming a gas with an adi- 
abatic index 5/3, the shock velocity with respect to the unper- 
turbed gas is V s = 4/3Vto P , whereas the pressure in the post- 
shock region (i.e., in the cocoon) can be expressed in terms 
of the velocity of the piston V p = V top and the upstream gas 
density po as (see, e.g., Kulsrud 2005) 



Jtop = ^ PoV s = - poV top . 



(9) 



By comparing this result with the condition of pressure bal- 
ance f top ~ Bq/Stt at the contact discontinuity at the top of the 
tower, we immediately see that the tower grows with a veloc- 
ity of order the Alfven speed computed with the unperturbed 
density po: 



top 



Va 



Bo 



(10) 



On the other hand, as we have shown earlier, y top should be 
of the order of the differential rotation speed /? Af2. Thus, 
we obtain the following estimate for the radius of the tower in 
terms of VPo, AS1, and po: 



Ro 



{ An 



1/3 



(47rpo) 



-1/6 



(ID 



We can also relate the radius of the tower to the radius of the 
central accretion disk. The poloidal flux can roughly be 
estimated as 

Vo~B d R d , (12) 



where B c [ is the typical magnetic field in the disk (or, rather, 
its dipole-like part), and Rd is the characteristic radius of the 
inner part of the disk — the base of the tower. Then, 
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Bo 



R d 



and 



Ro 
Rd 



Va.6 

v d 



1/3 



(13) 



(14) 



where Vd = RdAil is the characteristic differential rotation 
velocity in the disk (on the order of the Keplerian velocity Vk) 
and 

V A4 =-^L= (15) 
V'+Kpo 

is a convenient composite quantity that has a form of an 
Alfven speed involving the disk magnetic field and the unper- 
turbed star's plasma density; it doesn't represent a physical 
velocity and so can be arbitrarily high; in particular, it can be 
greater than the speed of light. 

Notice that as the tower makes its way through the star, 
two of the three main parameters, ^>o, and AO, remain un- 
changed, but the third parameter, the unperturbed stellar den- 
sity po measured at the top Z top of the tower, changes. In fact, 
it drops rather rapidly for a typical collapsar progenitor. Cor- 
respondingly, the radius of the tower increases as the tower 
grows. 

Also notice that, by assumption, the pressure spreads 
rapidly along the cocoon. Therefore, the characteristic mag- 
netic field and hence the radius of the tower are more or less 
constant along the tower at any given time. This justifies the 
assumption that the tower is a straight cylinder. In a more gen- 
eral situation, when the pressure equilibration in the cocoon 
lags somewhat, the radius of the tower will be a function of the 
vertical coordinate z. However, as is seen from equation (II 1L 

Ro scales only weakly with po (as p ^ 6 ). This fact enables us 
to regard the constant density case as a good approximation. 

2.4. Let's Do the Numbers 

Now we will use the above relationships to make some 
quantitative estimates based on the scaling relationships ob- 
tained in the previous section; thus, unavoidably, these esti- 
mates can only provide order-of-magnitude estimates. 

We assume that the core of the star has collapsed into a 
black hole of fiducial mass M = 3M with a corresponding 
gravitational radius R g = GMj c 2 ~ 5 km, and that some of the 
continuously infalling material has formed an accretion disk 
around the black hole. The pre-existing stellar magnetic field 
has been greatly amplified by the collapse because of flux- 
freezing and by the turbulent dynamo in the disk. Akiyama 
& Wheeler (2003) have argued that the field may reach the 
level set by the equipartition with the MRI-driven turbulence, 
as strong as 10 7 G. This field is mostly toroidal, how- 
ever. The large-scale poloidal magnetic field that we need in 
our model may require a large-scale helical dynamo (Black- 
man et al. 2006) for its production and will probably be some- 
what smaller than the toroidal field. Thus, we believe it is not 
unreasonable for the poloidal field at the disk surface to be a 
more modest B d ~ 10 15 G (see also Proga et al. 2003). 
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Now, if we take fiducial disk radius of Rj ~ R\sco(a = 
0) = 6R g ~ 3 • 10 6 cm, where a is the normalized Kerr pa- 
rameter of the black hole, we may have an initial poloidal 
dipole flux per unit toroidal angle of the order of = B^Bo ~ 
10 28 Z?i 5 /^ 65 in cgs units. (Note that if the black hole is 
spinning rapidly, then the inner disk radius is actually much 
closer, R\sco — R g — 5 km.) This poloidal flux is being con- 
tinuously twisted by the differential Keplerian rotation of the 
disk, with characteristic angular velocity Ail = il^Rd) — 

4-10 3 sec- 1 (M/3M Sun ) 1 / 2 ^ 3 6 /2 5 . 

The typical background density of the stellar material into 
which the tower propagates is po = 10 6 g/cm 3 . Hence, accord- 
ing to equation (I14> . the tower radius can be estimated as 



Bo ~ Rd 



V d 



= 3R d 



Bd,i 



Rd.b.sAili 



1/3 



(16) 



resulting in Bo ~ 10 7 cm for our fiducial parameter values. 

Using this estimate for the outer radius of the tower, we can 
get the following expressions for all the other parameters: 



V Afi = 



*0 

Rl 


■(£)" 




o.ifl d ( 


Rd, 6.5 AQ 


B±15 


Bo 


= 3- 10 10 





3.5 



2/3 



1/3 
P0.6 



10 14 G; 



(17) 



,1/3 d 2 / 3 AO 2 / 3 n-V^ID^ 
} dA5 K d,6.5 AS2 3.5 Po,6 



Notice that our crude estimate results in V top ~ Va,o being 
comparable to the speed of light c; moreover, as we shall show 
in § 0] a growth velocity of V top ~ 2.5Va,o is actually more 
typical. These facts strongly suggest that a fully-relativ istic 
treatment of the problem would be more appropriate (see § 15.11 
for discussion). Such a treatment lies beyond the scope of 
the current paper, but we do intend to develop the relativistic 
magnetic tower theory in a subsequent paper. 

Also, we can estimate the post-shock pressure in the hot 
cocoon above the tower as 



B, 



top 



A- 10 26 erg cm" 3 fl 2 



0.14 ■ 



(19) 



where Bo, 14 = Bo/(10 14 G). At such very high energy densi- 
ties the radiation pressure most likely dominates over the gas 
pressure; we can therefore estimate the plasma temperature in 
the post-shock region as 



'top 



( 3-Ptop 



a 
-15, 



1/4 



2-10 10 K~2MeV, 



(20) 



where a ~ 7.6 • 10~ l;> ergcm" 3 K" 4 . On the other hand, since 
we are dealing with a strong hydrodynamic shock between the 
cocoon and the unperturbed stellar material, the baryon den- 
sity in the cocoon is simply 4po, and so the baryon rest-mass 
energy density is 4,00c 2 ~ 4 • 10 27 erg cm" 3 po.6, and hence still 
exceeds the radiation/pair energy density by an order of mag- 
nitude. 

The total magnetic energy contained in the tower of 
height Z top can be estimated as 



s (0: 



■ 2irRoZ t0 p 



7T 



10 50 ergB 2 - 



>Bn 



(21) 



which represents a noticeable fraction of a typical GRB en- 
ergy (the factor 2 is added to take into account the fact that we 
actually have two towers, one above and one below the disk 
midplane). 



3. MATHEMATICAL FORMALISM 

The estimates presented in the previous section are very 
rough. A more rigorous mathematical model is needed to 
demonstrate how magnetic towers work and to obtain more 
accurate quantitative estimates of the tower growth parame- 
ters. We shall proceed to construct such a model. In partic- 
ular, we shall develop a general analysis in this section and 
then use it in §|4]to illustrate our model by specific examples. 

Consider a magnetic tower making its way through a star. 
Since the system is axisymmetric, the magnetic field can be 
written in cylindrical polar coordinates (R, <j>,z) as 



1 



pol 



B(B,z) = B pol + B^=-[Wx0] + 

K K 



(22) 



where ^(R,z) is the poloidal magnetic flux function, 7 po i(B,z) 

is (2/c) times the poloidal electric current, and <f> is the 
toroidal unit vector. If the magnetic field inside the tower is 
force-free, 7 po i is constant along poloidal field lines, I po \(R,z) = 

/polW- 

We shall focus on the middle part of Lynden-Bell's mag- 
netic tower (his section II), where the poloidal magnetic field 
is nearly straight and vertical. For simplicity we shall as- 
sume that the outer radius of the tower, Bo, is independent 
of height z. This means that, at any given moment of time t, 
all magnetic quantities depend only on one coordinate — the 
cylindrical radius B. Then, the force-free Grad-Shafranov 
equation — the main partial differential equation (PDE) that 
governs the field structure — reduces to a 2nd-order ordinary 
differential equation (ODE) for the function 'I'(B): 

' ' ' -/nolO*), 



Rd R 



(23) 



supplemented by two boundary conditions, 

\&(R = 0) = * = *(B = B ) . (24) 

As a side remark, we could consider a more general sit- 
uation where the outer radius slowly changes with height, 
Bo = Bo(z). Then, instead of just one ODE, we would have a 
one-parameter infinite set of ODEs, one for each height z. All 
these equations would have the same form as equation Q3t : 
the local height z would come in indirectly through the outer 
boundary condition ^> = ^o atR = Ro(z). 

Equation i23\ involves a free function, the poloidal cur- 
rent 7 po i( , i'). In reality, this function is not arbitrary but is 
determined by the differential rotation of the footpoints. By 
a simple geometrical consideration, the twist angle A$(5') 
(i.e., 2ir times the number of rotations) can be written as 



A$(*) = / po iW 



dlpoi 
R 2 B, 



pol 



: / P 0lW 



dz 



(25) 



where the integration is performed along the field line and / po i 
is the path length along the poloidal field. Assuming that most 
of the twist occurs in the main part of the magnetic tower, 
where the poloidal field is nearly vertical, we can write: 



A$(*) = / po i(*)Z top (0 



1 



1 



B 2 (*)|B zl (*)| B 2 (*)|B, 2 (*)| 



where Zt op (f) is the total height of the tower; its 



(26) 
time- 



dependence reflects the fact that the tower grows with time. 
Here, Rii^S) and R2(^) are the cylindrical radii of the two 
segments of the flux surface ^, and B z i ('J) and B^i^) are the 
values of the vertical magnetic field on these segments. 
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On the other hand, the twist angle is determined by the disk 
rotation law, 

A$(#) = Afi(10? . (27) 

where Af2( , J') is the difference between the angular veloc- 
ities of the two disk footpoints of the field line ^. Thus, 
determining the twist angle really requires the knowledge of 
both the radial disk rotation profile, fld(r), and the radial dis- 
tribution of the poloidal magnetic flux on the surface of the 
disk, ^d(f)- Whereas there exists a very natural choice for 
the first of these two functions — the Keplerian rotation law, 
fld(r) = f2#(r) ~ r~ 3 / 2 , the second function, i.e, ^d(r), is not 
very well known. In reality it is determined by complicated 
MHD processes of turbulent dynamo and turbulent magnetic 
flux transport, which are presently understood only poorly and 
definitely lie beyond the scope of the this paper. Therefore, 
instead of pretending that we know how ^/(r) is really deter- 
mined, we shall regard this function as essentially arbitrary. 
Correspondingly, instead of specifying it explicitly, we shall 
choose an indirect prescription dictated mostly by our mathe- 
matical convenience. In particular, we shall take ^(r) to be 
such as to result in a convenient and simple functional form 
for /pd ('J'). Thus, we shall use the following scheme: first, we 
shall pick a function /poK^P) that will make our calculations 
easier; then, we shall solve the Grad-Shafranov equation for 
this choice of Zpoi(^); and finally, we shall use equation i26\ 
to find ASl('I') and hence ^/(r) a posteriori. This approach is 
similar to what is known as the "generating-function method" 
in the theory of force-free equilibria (e.g., Uzdensky 2002). 
Even though the main logic of such a scheme is physically 
backwards, it makes sense to use it as long as we don't get 
an unreasonable functional form for ^j(r) in the end. At any 
rate, the alternative direct method, based on an explicit speci- 
fication of ^(r) and on a subsequent use of equation i26i to 
determine the corresponding / po i( l I'), would be more difficult 
to implement and, at the same time, would also lack a good 
physical justification since the prescription for ^(r) would 
necessarily be quite arbitrary. 

That said, there are some important features that this func- 
tion is expected to have. For example, we expect the separa- 
trix ^f s to correspond to the polarity inversion line on the disk 
surface, i.e., a line where B z d goes through zero. In general, 
however, the first radial derivative of B z (r) is not zero, and so 
generically, we expect the disk poloidal flux function to be- 
have as ^d — ^s + 0[(r-r s ) 2 ] near this point. Since the first 
radial derivative of flxif) is finite at this point, we may expect 



as * — > 



(28) 



Next, we need to determine the radius Rq of the tower. This 
is done by equilibrating the inside magnetic pressure to the 
outside cocoon gas pressure at Rq: 



Bl(Rp) | gfCgo) 
8tt 8tt 



side ■ 



(29) 



For simplicity we shall assume that the pressure is more 
or less uniform throughout the side part of the cocoon and 
constitutes a certain fixed fraction rj of the gas pressure in 
the post-shock region directly above the tower: Aide = ?7^top, 
77 < 1. The latter, in turn, is determined from the condition 
of force balance across the contact discontinuity at the top of 
the magnetic tower. Following Lynden-Bell (1996), we ig- 
nore the complications resulting from the complex shape of 
this boundary and treat the tower top as a solid circular lid 



of radius Rq. Then we only need to satisfy the integral force 
balance, which can be written as 



Ptop 



I 



(30) 



where F, is the total magnetic force acting on the top of the 
tower. It can be computed by integrating the vertical flux of 
momentum associated with Maxwell's stress over the tower's 
cross-section: 

ft] 



F 7 = 2tt 



8- 



8tt 



RdR. 



Using equation 1221 . we can rewrite this as 



ft> 



F, = ■ 



'pof 



dR 
~R~ 



(31) 



(32) 



Finally, once the pressure at the top of the tower is known, 
we can use equation (|9jl to estimate the speed at which the 
tower expands vertically, 



13 P, 



top 



(33) 



3.1. Formulation of the Problem in Dimensionless Variables 

In the following, it will be convenient to rescale R by 
the tower radius Rq, poloidal flux \fr by the total flux in the 
tower ^q, and poloidal current 7 po i by \&o /Ro- Thus we define 
new dimensionless variables 

A> (34) 



*0 



(35) 



(36) 



It will also be useful to introduce a characteristic magnetic 
field strength 

Bo = J; ■ (37) 

The Grad-Shafranov equation d23l > can be written in these 
dimensionless variables as 



Xl),.[ 1^ 



■n'm. 



(38) 



Geometrically, the magnetic tower consists of two regions 
(see Fig.|3}: 

1) inner region I of descending poloidal flux: < x < x s ; 
and 

2) outer region II of ascending poloidal flux: x s < x < 1. 
The boundary between these two regions is the separatrix 

surface x = x s . By assumption, each field line leaving the disk 
and rising up in region I will come back down to the disk 
through region II; therefore the total amount of poloidal flux 
in the two regions is the same. We shall need to solve the 
Grad-Shafranov equation J38l > separately in each of these two 
regions and then match the solutions at the separatrix. Corre- 
spondingly, we set the boundary conditions as: 

iP'( x = 0) = ^"( x = l)= l, (39) 

^{x = x s ) = ^"{x = x s ) = ^ s . (40) 
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FIG. 3. — A schematic drawing of a magnetic tower. The lines attached 
to the disk represent poloidal magnetic field lines (toroidal magnetic field is 
also present). 



Without losing generality, from now on we shall set ip s = 0. 

The position of the separatrix x s is determined from the 
force-balance condition across it: 



(B^+B$' = (B^+Bp" (41) 

This condition can be greatly simplified, however. Indeed, 
first, because = I po i(ip)/R and because the two sides of 
the separatrix correspond to the same field line ip = ip s , and 
hence l'(ip s ) = l"(il>s), we have (B 2 ^)' = (B 2 ^)". (Moreover, we 
shall actually assume that / — > ^/tp — > near the separatrix, so 
B,/) actually goes to zero at the separatrix.) Thus, our force- 
balance condition becomes: 



(B 2 z y = (B 2 ) 11 = 

The vertical force F z is 
F z = 4 B o R o 



[/ 2 (V0- 



at x = x t 



dx 



and hence the pressure at the top of the tower is 

l 



top ■ 



1 9 

Br 

4vr c 



2 dx 



(42) 



(43) 



(44) 



The condition ( I29i of pressure balance across the side wall 
of the tower, x = 1, then becomes 



I 2 Oj = 1) + ip](x = 1) = 8tt = 8vrr/ 

B n 



1 top 



1 



--2 V J[I 2 W)-ip 2 ]^, (45) 
o 



Finally, we can rewrite the expression i25\ for the field-line 
twist angle A <£>('!') in our dimensionless variables as 



V, 



1 



1 



(46) 



_xi\if; x (xi)\ x 1 \i\) x {xi)\_ 

where x\ = x\ (ip) and X2 = X2(ip) are the radial positions of the 
inner and outer segments of field line ip, respectively. 

4. EXAMPLES 

When dealing with any complicated physical system, it is 
often useful to have at hand a few analytical examples that 
could be used to explicitly illustrate the general characteristic 
behavior of the system. In our particular problem, the simplest 
way to obtain explicit analytical solutions is to consider the 
cases in which the Grad-Shafranov equation fl38l > becomes 
linear. The most general functional form of /(i/;) that leads to 
a linear Grad-Shafranov equation is 



I(tp) = y/aip 2 + bijj+c. (47) 

For simplicity, we shall consider the situation where the 
poloidal current vanishes along the separatrix, = 0) = 0, 
which makes c = 0. Other choices for c are possible and may 
result in interesting solutions. For now we consider c = and 
thus have _ 

I(ip)=^VaJTb. (48) 

There are several interesting special cases that we will con- 
sider in this section. 

As we shall see, these examples should be seen as math- 
ematical idealizations. In particular, the specific functional 
form of I(ip), given by equation J48b . may lead to solutions 
that are not regular at the axis. A physically-motivated regu- 
larity requirement at R = means that I(t/j) should go to zero 
at the axis fast enough (e.g., like 1 - ip). 

Mathematically, solutions obtained in this section can be 
seen as describing the magnetic field behavior well outside 
of a narrow axial boundary layer, the inner "core jet", inside 
of which the poloidal current distribution deviates from 
equation J48b . The resulting singularity in our Case 2 in par- 
ticular is similar in nature to the divergent behavior of the az- 
imuthal magnetic field of a line current carried by a thin wire 
(inside of which the magnetic field is of course regular). 

4.1. Case 1: a Configuration without a Central Line Current 
First, we shall consider the case in which I(ip) behaves 



as near the separatrix -0 = and at the same time goes 
to zero at the rotation axis: 7(1) = 0. Thus we choose b = -a = 
2 > in equation (I48> and hence 



which corresponds to = /i 2 (ip-l /2). 



(49) 



By substitutions tp = xu(x) + 1 /2, this equation is reduced 



to the Modified Bessel's equation x u +xu - u(fi x + 1) = 0, 
and so the general solutions in the two regions can be written 
in terms of modified Bessel's functions as: 

">pi{x)= - + a\xI\(px) + b\xK\(p,x) , (50) 



2 

^>n{x) = — + a2xl\ (px) + b2xK\ (fix) . 



(51) 



Here, ct\, b\, a2, and Z?2 are the arbitrary integration con- 
stants. They are determined by the boundary conditions j39l>- 
< I40I that can be cast as a linear system of four equations. Solv- 
ing this system we get the explicit expressions for these coef- 
ficients in terms of the parameter /i and the position of the 
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separatrix*,: 



a\ ■ 



bi = - 



i+yKdy) 
2yh(y) 



1 

a2 = -7T 



b 7 = — 



yKi(y) + fiKi(n) 

2y /i(y)Xi(Ai)-Zi(y)/i(A*) 
1 yh(y) + jih(ji) 



2y h<y)KiQi)-Ki(yyii(jMy 



(52) 
(53) 
(54) 
(55) 



where we defined y = jjx s to simplify notation. 

Next, the separatrix position x s as a function of u is deter- 
mined by the separatrix force-balance condition I42i that can 
be written as 



(a i + a 2 )l[{y) + (b, + b 2 )K[{y)= ^ 



r 



(56) 



Thus, we get a non-trivial transcendental algebraic equation 
that determines y as a function of ji. Once the solution of this 
equation is found, we can use equations (1551 to calculate the 
entire magnetic field structure. We solved equation d56t nu- 
merically using Mathematica. We found that the resulting de- 
pendence y(p) is very close to linear y = fi/y/2 for small and 
finite values of \i. For ji ^> 1 the slope of the linear depen- 
dence slightly changes and the asymptotic behavior becomes 
yip — > oo) ~ fj,— 7. Correspondingly, the function x s {p) varies 
from 1 / y/2 for /_t — > to 1 for \i — > oo. [The small-/i behav- 
ior is very easy to understand: it corresponds to the situation 
in which the toroidal field is negligibly small, and so the hori- 
zontal force balance dictates that the magnitude of the vertical 
field be uniform inside the tower. Since the total upward verti- 
cal flux must be equal to the total downward flux, this means 
that the area of region I must be equal to that of region II; 
hence the separatrix between them has to be at x s = 2 -1 / 2 .] 

Finally, we need to fix the value of \i. We do this by us- 
ing the condition J29i of pressure balance between the mag- 
netic field in the tower and gas pressure in the cocoon that 
surrounds and confines the tower. As discussed in § |5J this 
condition can be viewed as the condition that determines the 
tower's radius once the poloidal current 7 po i is specified. Here 
we, in effect, turn this around and determine fi (which sets the 
overall scale for 7 po i) in terms of the tower's radius. 

The solution of course depends on the assumed value of rj. 
For a given rj, one can obtain a series of solutions for various 
values of /i and then find the value no(rj) for which condi- 
tion J29l > is satisfied. We have done this using Mathematica 
and found that the allowable range of rj is limited. Namely, 
solutions exist only for r\ > <q c w 0.08, with jJ,o(rf) diverging 
as {rj-rj c )~ l in the limit rj — > rj c . Both P top and V top also di- 
verge and x s approaches unity as 77 — > ij c . 

We then investigated a few values of r\ in more detail. For 
77 = 1 we found (1q(j] = 1) ~ 6.70 and x s (r] = 1) ~ 0.705 (just 
a little less than 1 / \2). The corresponding magnetic field 
components as functions of x are plotted in Figure|4](the top 
two panels). We see that both the vertical and toroidal field 
components increase towards the axis, with ~ x~ l l 2 . At 
the separatrix (marked on the plots by the vertical dashed 
line), the vertical field, as expected, reverses, whereas the 
toroidal field goes to zero on both sides of the separatrix as 
|jc— Xs] 1 ' 2 . The post-shock pressure above the tower is found 
to be P t0 p(i] = 1) — 0.945 fig and the vertical expansion speed 
of the tower, given by equation i33\ . is V top (rj = 1) ~ 3 Va,o- 





Fig. 4. — Vertical (left) and toroidal (right) magnetic field components for 
two solutions in Case 1 , computed for r) = 1 .0 (top panels) and rj = 0.3 (bottom 
panels). The vertical dashed lines on each plot shows the position x s of the 
separatrix, across which the vertical field component (B z ) reverses sharply 
[the toroidal (Ba) component is symmetric with respect to jc s ]. 



It is interesting to note that, as can be seen in the top 
left panel of Figure |4] the vertical magnetic flux tends to 
be pushed out (by the toroidal field pressure) from the 
intermediate-radii parts of each of the two tower regions 
(x ~ 0.4-0.5 in region I and x ~ 0.8 in region II). As rj de- 
creases, this effect becomes even more pronounced, as can 
be seen from the bottom left panel of Figure |4] which corre- 
sponds to rj = 0.3. Also, in agreement with an earlier dis- 
cussion, the pressure at the top of the tower, ftop, expressed 
in units of Bq, and the tower growth velocity, V top , expressed 
in units of the Alfven speed Va.q, both grow rapidly as 77 de- 
creases. For example, in the case rj = 0.3, presented in the 
bottom two panels of Figure^ we find P to p(?7 = 0.3) ~ 8.4B, 2 
andV top (77 = 0.3)~8.9V A , - 

4.2. Case 2: a Configuration with a Central Line Current 
Let us now consider the case a = 0. Then, 



1(f) = ^bip ee \y/2i{; 



(57) 



Notice that the poloidal current in this model does not van- 
ish at the tower's rotation axis ip = 1, implying the existence 
of a non-zero axial line current. Such a current may repre- 
sent, for example, a very narrow core jet surrounded by the 
magnetic tower, perhaps produced by the Blandford-Znajek 
(1977) mechanism acting along the field lines that thread the 
black hole. This line current may have important implications 
for the tower's growth speed and collimation, as we will dis- 
cuss below. 

The Grad-Shafranov equation J38i becomes: 



xi), [ — ) = -A" = const . 



and can be easily integrated: 

„ - 2 

ip(x) 



\ lX - logx+Ci-+C 2 . 



(58) 



(59) 



This expression represents the general solution of the 2nd- 
order differential equation (1581 . It involves two arbitrary in- 
tegration constants, C\ and C2, that are determined by the 
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FIG. 5. — The function A(,r s ) corresponding to equation 1621 for our Case 2. 

boundary conditions J39I — d40i . The resulting final expres- 
sions for ijj(x) in the two regions are: 

^( x ) = -^-x 2 \ g(-) + l-^, (60) 

L \X S J X s 

X 2 ( 1-x 2 \ x 2 -x 2 

ip"(x) = - — f x 2 logx - — 2 x] logx s ) + -j— f .(61) 

Finally, the position x s of the separatrix is obtained by 
substituting these solutions into the separatrix force-balance 
condition J42i . yielding a transcendental algebraic equation 
forx s (A). This equation can actually be resolved with respect 
to A as a function of x s , resulting in an explicit expression: 

-t/rfeP-- (62) 
x s y i-xf+iogx s 

This function is plotted in Figure[5] One can see that there are 
two separate allowed ranges of x s : < x s < x\, and X2 < x s < 
1, where w 0.45076, and x 2 = 1 /a/2 w 0.707. At x min « 
0.278 the function has a minimum A(x m ; n ) w 7.824. Forx, > 
X2, the solution X(x s ) increases monotonically starting from 
A = at x s = X2 and diverges as A ~ [2/( 1 - x s )] 1 1 2 as x, — > 1 . 

The next step in our program is to use the obtained solution 
to determine the vertical expansion velocity of the tower, V top . 
To do this, we first need to compute the integrated magnetic 
stress on the tower's top, given by equation d43i . Notice, how- 
ever, that in this model, as in any model with a non-zero line 
current along the axis, the toroidal magnetic field diverges as 
1/x as x — > 0, and hence its pressure diverges as 1 /x 2 . Then, 
the vertical force due to the toroidal field pressure integrated 
over the top lid diverges logarithmically at small radii. To 
manage this singularity, we will modify the problem some- 
what by introducing a small inner cut-off radius R m = x m Ro of 
the tower, x; n <C 1 . This inner radius corresponds to the radius 
of the central jet that is embedded in the tower and that is ac- 
tually responsible for carrying the axial line current. We use 
Xj n as a small parameter to facilitate our analysis. 

The contribution of the toroidal field pressure to the overall 
integrated vertical force, 

F^^RlHxdx=^ 2 [l\i,A, (63) 
J 8vr 4#g J x 

is dominated by the contribution near the cut-off radius; using 
the function given by equation d57i . this contribution can 
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FIG. 6. — The position x s of the separatrix as a function of r] for the four 
solution branches corresponding to .Vj„ = 10T 4 in our Case 2. 

be estimated as 

l 

= Y§J ^dlogx ~ | lB 2 R 2 . (64) 

An 

where we introduced a large parameter 

7= |logx in |> 1. (65) 

Now, for what values of the parameters rj and 7 do solutions 
exist? To answer this question, we have used Mathematica 
to map out the two-dimensional (77,7) parameter space. We 
have restricted our investigation to 7 ^> 1 and rj < 1. We 
have found that the system exhibits a very rich mathemat- 
ical behavior characterized by the multiplicity of solutions. 
The overall picture is illustrated in Figure [6] where we plot 
the separatrix position x s as a function of rj for a fixed value 
of 7 = 9.21 (corresponding to x; n = 10~ 4 ). We found that es- 
sentially there are two pairs of solutions, one pair correspond- 
ing to x s < X2 = 2" 1 / 2 and the other corresponding to x s > Xz- 
When 77 = 1 there are only two solutions, one in each of these 
two regions. However, for 77 < 1 two more solutions appear, 
so there are four in total, with very different properties. As 
one decreases 77 at a fixed 7, the difference between the two 
solutions in each pair gradually diminishes, until they finally 
merge into one solution. As can be seen from Figure [6] this 
happens at two critical values of 77, which we call 77mi n ,i for 
the two solutions below X2 and ?7 m j n ,2 for the two solutions 
above Xz- We found that both ?7 m i n ,i and ?7 m i n .2 scale inversely 
with 7: 77 min j(7)w2.257" 1 ,and?7 m i n ,2(7)« I.57" 1 . There are 
no solutions for 77 < ?7min, 1,2(7) • 

As an illustration, we present the four solutions correspond- 
ing to x; n = 10 -4 (7 = 9.21) and 77 = 0.3 obtained using Math- 
ematica. In Figure we plot the vertical and toroidal com- 
ponents of the magnetic field for these solutions. One can 
see that the toroidal field goes to zero as |x-x s | 1//2 on both 
sides of the separatrix without changing sign and grows as x -1 
near the axis, which is explained by the assumed shape of 
the function I(tp). One can also notice that the vertical field in 
these solutions not only reverses sharply across the separatrix, 
but also changes sign by going smoothly through zero at one 
or two other points. This means that the magnetic field de- 
scribed by these solutions actually has a more complex topol- 
ogy than that of the simple tower of Figure [3] Thus, solu- 



12 



b) 




d) 



FIG. 7. — Vertical (left) and toroidal (right) magnetic field components of 
the four solutions in Case 2, computed for x m = 10 -4 (7 ~ 9.2) and r) = 0.3. 
The vertical dashed lines on each plot shows the position x s of the separatrix, 
across which the vertical field component (B-) reverses sharply [the toroidal 
(B^) component is symmetric with respect to x s ], 



tions 2c and 2d correspond to two towers, one inside the other, 
and solutions 2a and 2b corresponds to three nested towers. 

The reason for the existence of two solutions with x s > X2 
for 77 > 77 m i n ,2 can be traced to the following argument. The 
vertical force J64i leads to a cocoon side pressure that scales 
as A 2 7?7Bo- This pressure has to be balanced at the outer edge 
of the tower by the internal magnetic pressure that generally 
is of the order of Bq/(1 -x s ) 2 . It turns out that, provided that 
V > f?min,2, there are two ways to make such a balance, two 
different regimes, distinguished by whether x s is close to 1 or 
not. We shall now consider these two possibilities (which we 
call Case 2c and Case 2d) separately and give basic analytical 
derivations under the assumption that 777 3> 1 (and hence 77 3> 

?7min,2). 

Case 2c: Let us assume that x s — 1 is finite. Then, the force 
balance condition at the side wall of the tower can be satis- 
fied only if A is small, namely, A 2 oc (777)- 1 . This suggests the 
following overall picture. Since A 2 is the factor that parame- 
terizes the general strength of the toroidal field pressure, we 
see that the toroidal field term — 77'(\&) in the Grad-Shafranov 
equation is unimportant in most of the tower, with the excep- 
tion of a small vicinity of the inner edge. The field is then 
predominantly vertical almost everywhere; the force-free bal- 
ance then dictates that this field must be essentially uniform 
in each of the two regions. Since the total upward-directed 
magnetic flux in region II must equal to the total downward- 
directed flux in region I, the areas of the two regions must then 
be the same, hence x 2 « 1 /2. Then we can approximate x s in 
equations ( I60t -( I6H by x% = 2" 1 / 2 and rewrite them as 



i/j'(x)~1-2x 2 



(66) 



7/> /7 (x)-2x 2 -l : 



(67) 



corresponding to Bi(x) ~ -4Bq, and B''{x) ~ 4Bq. 

Next, the overall integrated force F z is the sum of the pos- 
itive contribution of the toroidal field pressure near the 
axis and the negative contribution from the poloidal field ten- 
sion in the rest of the tower: 



F=F 7 



1 



■ B l R l 



ft 



dx 



A*7 
2 



BnRn . 



(68) 



Then, using the condition ( 145 \ of the pressure balance be- 
tween the tower and the cocoon at the outer edge x = 1, we 
obtain a simple approximate expression for A in terms of 7: 



A~2, 



I + 77 

777 



< 1. 



(69) 



which reduces to A ~ (8/7) 1 / 2 for 77 = 1 and A ~ 2(r\ r ff x l 1 
for 77 -c 1 . By combining this expression with the asymptotic 
behavior 



\(x s —>x 2 )~4-\ 



V2 



1-W2 



\/V2~%.5^{x s -\/V2) 1 



n 



(70) 

of the function \(x s ) [obtained from equation \62l in the limit 
x s — > X2], we can determine exactly how close x s has to be 
to X2 in terms of 7 and 77: 



-X 2 = ■ 



I+77 l-log2 
VI 4V2 



: 0.054 



l+T] 

?77 



(71) 



This becomes x s - X2 
-1 



0.1 17 for 77 = 1 and x s — x% ~ 

0.054(777)-' for 77 < 1. 

Finally, by substituting equation ( I69> into d68l > and using 
relationship J30t . we get an estimate for the cocoon pressure 
at the top of the tower: 



top 



-Bi 



0.637 



-B 



7T7] 



n 



• 



(72) 



which is of course consistent with the pressure balance across 
the side wall between the cocoon and the magnetic field 4Z?o 
inside the tower. Equation J33l > can then be used to estimate 
the tower's vertical expansion velocity: 



:2.45t7 -1/2 V, 



A,0 ■ 



(73) 



This result, together with its counterpart from the previous 
example, indicates that V top between 2.5 and 3 Va,q provides 
a more accurate estimate for the tower growth velocity than 
equation dlOt . 

With these results at hand, we can go back and compute 
the differential rotation profile ACl(ip) that corresponds to this 
case. As l ong as we are not too close to ip = 1, we can use 
equations i66\ and J67t to explicitly express x\(ip) and X2(i[f) 
in terms of ip and we can also approximate ip' x = -4x,jpU = Ax. 
By substituting these approximations, along with d57i . into 
equation J46l >. we get a very simple relationship 

t Rq 1 - Ip 1 



(74) 
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This expression breaks down near ip = 1 as of course it should 
since the toroidal field's pressure becomes significant and 
modifies the poloidal field structure near the inner edge. Also 
note that the factor V2X ~ 4| logjCi n |~ 1//2 that is in front of this 
expression, although formally vanishing in the limit x m <C 1, 
in reality is never very small; for example, it is equal tow 1.1 
forjc in = 1 x 10" 6 . 

We have used Mathematica to obtain exact (i.e., without 
relying on the assumption 7 ^> 1) solutions for this case for 
several values of 77 and 7. In particular, we found that, in the 
limit of large 7, the resulting 7-dependences of A, the separa- 
trix position x s , and the pressure in the cocoon P top at a fixed 77 
are in a good agreement with the above simple estimates. In 
addition, Figure shows the radial profiles of the vertical 
and toroidal magnetic field components inside the tower for a 
representative case Xm = 0.0001 and 77 = 0.3. One can see that 
the vertical field is indeed close to a constant 4Z?o outside the 
separatrix and -4Z?o inside, as expected, although it deviates 
significantly from this constant close to the central axis. 

Case 2d: Another interesting case is that corresponding 
to x s being close to 1 and hence a large A. Then, assuming that 
the contribution J64i of the toroidal pressure to the integrated 
vertical force dominates over the poloidal field contribution 
(this assumption is valid as long as 77 -C 1, as we shall show 
below), we can estimate the post-shock pressure at the tower 
top as 



Ptop — 



irR 2 0Ut 2tt 



(75) 



Applying the condition of pressure balance across R = Rq, 
we then get: 



1) * 2 n 



77 + i>l(* = 1) = 87TP sld e = STT^top = 4A 2 7 77^ . 

tin Rn 



(76) 

The contribution of the toroidal field to the pressure balance 
at the outer boundary is B 2 Jx = 1) = = = 2A 2 fi 2 >. 
Assuming that 777 ^> 1, we can therefore neglect it in compar- 
ison with 87rP s ide, and so the above pressure balance condition 
yields 

^=1)~2AVT7- (77) 

Now we want to use this relationship in conjunction with 
our general solution derived above to determine the param- 
eter A. Notice that the introduction of the inner radius does 
not alter the general functional form of ip(x) given by equa- 
tion ( I59> . The inner boundary condition is changed into 
ip(Xm) = 1. However, we expect that, for x- m <C 1, the relation- 
ships between the coefficients that describe the function ip(x) 
are changed only slightly, even though they had been derived 
using the boundary condition ip(x = 0) = 1 . Then, as long as 
we are not too close to the inner boundary, we can still use 
expressions 

d6QH6D for our solution. In particular, we can 
write tp x (x = 1 ) as 

Ji x 2 



ip x (x= 1) = 



1- 



■ A 



1- 



logXs- — - 



(78) 



We can simplify this expression by noticing that the rela- 
tionship J62l > between A and x s can be cast in the form 
A 2 x 2 logx, = 2(2x 2 - 1)- A 2 x 2 (l -x 2 ). Then, the above expres- 
sion for ip x (x =1) greatly simplifies and, substituting it into 
the pressure balance condition ( 1771 . we obtain 



tj, x (x= 1) = 4+A 2 



:2A 



777. 



(79) 



777- 



A 2 



(80) 



From this, we can express x s in terms of A as 
2 1 „ A 

This equation, together with equation ( I62> . completely deter- 
mines both x s and A for a given value of the product 77; 3> 1 . 
We can go even further and get explicit formulae for x s and A 
by making use of the fact that x s » 1 [and hence, according to 
eq. i62i , A 1] in the regime under consideration. Then, the 
first term in i79i is small compared with the second term and 
we immediately obtain 



A ~ 4V7?7 = 4v/?7|logjq n | > 1 , 



(81) 



I logx,| ^ 1 -x s ~ ^ ~ J- < 1 . (82) 

With these results, we can now go back to our solution 
given by equations ( I60> - (I6U and derive the following asymp- 
totic expressions for ip(x) in the two regions 



lf/(x) ~ 1 — 2x 2 — 8777 X 2 logX , X 3> Xi n ; 



(83) 



^"(x) ~1 -8777(1 -x 2 +x 2 \ogx)~ 1 -8777(1 -x). (84) 

From our estimations of the vertical force and the post- 
shock pressure above the tower, we find 



7A 2 



BlRl-SYvBiRi, 



2„ j}2 r,2 



87 2 , 



Correspondingly, the vertical expansion speed is 



(85) 



(86) 




2477 7 V a ,o»Va,o. (87) 



This suggests that magnetic tower configurations with a non- 
zero line current along the axis may be especially effective in 
producing jets that are rapidly propagating through a confin- 
ing external medium. 

As a final check, we now go back to our estimate of the 
integrated vertical force and ask under what circumstances it 
is indeed dominated by the toroidal field contribution. Using 
our solution (I83> - (I84> and assuming that 7?/ ^> 1, the poloidal 
field contribution can be estimated as F! B ) ~ -&-f 2 T] 2 B^Rq, 
whereas, substituting J8 It into equation J64i . one estimates 
the toroidal field contribution to be F^^ ~ 87 2 77fi 2 1 /? 2 >. Thus, 
we see that F? B ^ indeed dominates as long as 77 -C 1. 

5. DISCUSSION 

5.1. Transition to the Relativistic Expansion Regime and the 
Final Opening Angle of the Jet 

As we have shown in § 12.21 the radius Rq of the tower, 
and hence the rotational velocity of the field lines and the 
tower's growth velocity Vt op ~ Af2/?o scale with the back- 
ground density as p '^ 6 . Therefore, as the tower expands into 
the outer, less dense regions of the star, Vtop inevitably reaches 
the speed of light at a certain point that corresponds to some 
critical density p re i. Our non-relativistic model becomes be- 
comes invalid at this point and any further expansion of the 
tower requires a relativistic treatment. For our fiducial val- 
ues B d = 10 15 G, R d = 3 • 10 6 cm, and AO, = 3 ■ 10 3 sec" 1 , 
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this transition to the relativistic regime takes place at a crit- 
ical density po.rei ~ 10 6 g/cm 3 [see eq. ( I18H . For a typical 
massive stellar GRB progenitor this density corresponds to a 
distance Z le i from the center that is on the order of a few times 
10 8 cm; according to equation d!6i . the radius of the tower 
corresponding to this critical density is R re i ~ 10 7 cm. 

A relativistic generalization of the magnetic tower model 
is beyond the scope of this paper but is something we in- 
tend to develop in the near future. Before a complete so- 
lution is achieved, however, we cannot really say anything 
definite about further expansion of the tower. However, we 
can, perhaps, look at the results of fully-relativistic hydro- 
dynamic simulations (Zhang et al. 2003) for some physical 
insight. In those simulations the relativistic jet remained col- 
limated as it propagated through the star, across several or- 
ders of magnitude in stellar density. The authors of that pa- 
per have attributed the observed collimation to recollimation 
shocks in the cocoon and to relativistic beaming in the jet. But 
the physical processes in the cocoon should not change if we 
replace the inner relativistic hydrodynamic jet with a relativis- 
tic magnetic tower. Furthermore, we expect some additional, 
magnetic collimation due to the toroidal-field hoop stress (the 
pinch effect). The bottom line is that we expect the magnetic 
tower to remain well collimated even after it transitions into 
the relativistic regime. In particular, we suggest that the fi- 
nal opening angle of the magnetically-dominated outflow will 
be no larger than the inverse aspect ratio of the tower at the 
moment of relativistic transition: 



~ Zrel 

Whether this prediction is true will have to be determined by 
a fully-relativistic analysis and by relativistic MHD simula- 
tions, which we hope will be completed in the near future. In 
addition, we acknowledge that some first steps toward devel- 
oping an analytical understanding of magnetically-dominated 
outflows in the ultra-relativistic regime were taken recently by 
Lovelace & Romanova (2003) and by Lyutikov (2006). 

5.2. Tower Structure 

Figure @] shows the run of vertical and azimuthal magnetic 
field components for Case 1 as a function of scaled radius x. 
For the rj = 1.0 case (bottom panel), the toroidal field scales 
as x~ l l 2 so that the enclosed toroidal magnetic energy in- 
creases linearly with radius. The toroidal magnetic energy is 
not concentrated toward the axis, but is distributed throughout 
the tower. For Case 2 with a central line current (Figure [7J, 
however, the toroidal magnetic field increases toward the axis 
(x = 0) as x~ l . The magnetic energy density is thus strongly 
peaked toward the axis for this case. This demonstrates that 
the presence of a central line current results in a more tightly 
beamed outflow. 

For the field structures of Case 2 plotted in Figure [7] 
we found solutions corresponding to multiple nested towers. 
Case 2c, for example, has B ; passing through at x ~ 0. 1 and 
a separatrix at x ~ 0.75 corresponding to two distinct tow- 
ers, one within the other. Field lines in the inner tower are 
tied to the disk deeper in the potential of the black hole or 
neutron star where the characteristic velocity becomes com- 
parable to the speed of light. Thus, the inner tower may be a 
more energetic component of the outflow and may reach the 
ultra-relativistic expansion regime earlier. In addition, the in- 
ner tower is shielded from the cocoon material by the outer 
tower which envelopes it. This is beneficial for achieving rel- 
ativistic flow in the inner tower since baryon contamination is 




j x=1 

FIG. 8. — A schematic drawing of a twin magnetic tower. The outer tower 
which forms a hollow cylindrical shell which surrounds the inner tower. 



suppressed. We note that, in principle, the tops of the towers 
can propagate at different speeds, as indicated schematically 
by the tower heights in Figure [8] However, we have not an- 
alyzed this case in detail since all analysis we present here 
assumes a single tower height. 

The physical picture presented in this paper, with a smooth 
coherent magnetic structure of the outflow, is an idealization, 
necessary for obtaining a basic physical insight into the sys- 
tem's dynaimcs and to get the main ideas across in the clear- 
est way possible. The actual magnetic field, especially if it 
is produced by a turbulent dynamo in the disk, is likely to be 
different from such a simple system of nested axisymmetric 
flux surfaces. Instead, it may consist of an ensemble of loops 
of different sizes and orientations. It may thus have a highly- 
intermittent substructure on smaller scales, both temporal and 
spatial. In fact, our Case 2 solutions discussed above, with 
two magnetic towers, one enveloped by the other, represent a 
simple specific example of this general situation. 

What is important though, is that each of these smaller mag- 
netic structures is subject to the same physical processes as a 
single large tower: twisting due to differential rotation and a 
subsequent inflation, controlled by the external pressure of the 
cocoon and of other loops growing at the same time. As a re- 
sult, a single tower may be replaced by a train of spheromak- 
like plasmoids, pushing each other out along the axis. The 
hoop stress still works inside each of them, and so the overall 
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dynamical effect may be the same as that of a single tower, at 
least qualitatively. The same picture is expected to develop if 
a single large-scale tower becomes unstable to reconnection 
and brea ks up into smaller plasmoids, as we discuss in Sec- 
tion !5.5l In either case, the resulting multi-component struc- 
ture of the outflow may be responsible for the observed inter- 
mittancy in GRBs. 

5.3. Baryon Contamination from the Disk 

The tower plasma may be polluted with baryons from the 
disk by two main processes, magneto-centrifugal acceleration 
(Blandford & Payne 1982) and radiative neutrino-driven ab- 
lation (Qian & Woosley 1996; Levinson 2006). In the case of 
magneto-centrifugal acceleration, the field lines must be in- 
clined outward by at least 30 degrees from vertical for a disk 
wind to be launched. It is of interest to note that, for our 
configuration, field lines forming the inner tower are plausi- 
bly constrained to leave the disk at progressively steeper an- 
gles, moving inward toward the rotation axis. This field con- 
figuration would suppress the magneto-centrifugal launching 
mechanism for the inner tower. The outer tower may be more 
baryon loaded than the inner tower and would constitute a 
less relativistic magnetically dominated flow at larger angles 
from the pole. Yet farther from the pole, the cocoon represents 
an even higher baryon loaded component and is consequently 
expected to be slower and contain a relatively small magneti- 
zation (whatever magnetic field is mixed into the cocoon by 
plasma instabilities). Future relativistic magnetohydrodynam- 
ical simulations will be able to address these issues and in 
particular determine the angle of field lines emerging from 
the inner disk enabling better estimates of baryon contamina- 
tion. In addition, we note that in the case of black hole-disk 
magnetic coupling, one should expect no baryon loading at all 
along the part of a field line that connects to the black hole, as 
has also been pointed out by e.g. Levinson & Eichler (1993). 

Some baryons are expected to be driven into the tower 
(along field lines attaching to the disk) by neutrino ablation 
(e.g., Levinson 2006). As of yet, no detailed calculations of 
neutrino-driven ablation exist for collapsar disks. For exam- 
ple, neutrinos may escape the disk in hot magnetically dom- 
inated bubbles without advecting large numbers of baryons 
into the tower. In addition, Vlahakis et al. (2003) argued that 
neutrons, which may constitute the majority of the baryons 
in the disk wind, are not subject to electromagnetic accelera- 
tion and hence decouple from the relativistic outflow at rather 
moderate Lorentz factors; this significantly lessens the baryon 
pollution expected from the disk wind. 

5.4. Effect of MHD Instabilities on the Tower Evolution 

One of the most serious uncertainties in the magnetic tower 
model is the stability of this highly twisted structure. In this 
section we discuss the effects that ideal MHD instabilities 
may have on the magnetic tower in the collapsar context. A 
scenario involving reconnection inside the tower due to finite 
resistivity will be discussed in the next section. 

Indeed, as the tower grows, it becomes very elongated and 
hence prone to a number of ideal MHD instabilities. Among 
the instability candidates that first come to mind there are: 

(i) Rayleigh-Taylor (or, perhaps more relevant, Richtmyer- 
Meshkov) instability taking place at the tower's top. This may 
cause splitting of the magnetic tower into a number of sepa- 
rate strands interlaced with regions filled with stellar matter. 

(it) Kelvin-Helmholtz instability at the side boundary be- 
tween the tower and the cocoon, leading to entrainment of the 



baryonic stellar material into the tower. This instability how- 
ever should be stabilized by the strong vertical magnetic field 
in the tower. 

(Hi) A non-axisymmetric ideal-MHD kink-like instability 
taking place in the main body of the tower. Although this 
problem has not been studied in the magnetic-tower context, 
it has been considered in the framework of relativistic force- 
free electrodynamics for field lines connecting a disk to a 
black hole (Gruzinov 1999). In addition, an important clue 
can be derived from stability studies of regular MHD jets. 
For example, a recent numerical investigation by Nakamura & 
Meier (2004) indicated that, when a non-relativistic Poynting- 
flux dominated jet propagates through a stratified external 
medium, the jet stability strongly depends on the background 
density and pressure profiles along the jet. In particular, these 
authors found that a steep external pressure gradient forestalls 
the instability onset. 

If it does go unstable, the kink is probably the most dan- 
gerous instability that may lead, in the nonlinear phase, to a 
significant, although perhaps temporary, disruption. Such a 
disruption, however, is not necessarily a bad thing: the tower 
may be able to reform after being disrupted (as is seen in lab- 
oratory experiments by Lebedev et al. 2005) and the resulting 
non-steady evolution may provide a plausible mechanism for 
rapid variability seen in gamma-ray bursts. 

We also would like to point out that, as one expected out- 
come of the nonlinear development of an MHD instability, 
a significant fraction of the toroidal magnetic field energy 
may be dissipated into thermal energy (Eichler 1993; Begel- 
man 1998). Importantly, this process may occur somewhere 
inside the tower, way above the inner accretion disk. Since 
the plasma density in the tower is very low, this may have 
an additional beneficial side-effect on the explosion dynam- 
ics due to the possibility of non-thermal neutrino emission 
coming from the magnetic dissipation sites. The resulting 
high-energy neutrinos are expected to have better coupling to 
baryonic stellar matter, thus providing an extra source for the 
strong shock propagating through the star (Ramirez-Ruiz & 
Socrates 2005). In addition, as has been shown by Drenkhahn 
& Spruit (2002; see also Giannios & Spruit 2006), magnetic 
energy dissipation may also have a positive influence on the 
acceleration of the Poynting-flux dominated outflow. 

So far as we know, there have been no formal stability 
studies of magnetic towers to date. Such studies are clearly 
needed. They may involve a linear perturbation analysis or 
a non-axisymmetric MHD simulation. Such stability stud- 
ies should take into account several stabilizing effects. First, 
Tomimatsu et al. (2001) have conducted a linear stability 
analysis of the kink instability in a narrow rotating relativis- 
tic force-free jet and have shown that rapid field-line rotation 
inhibits the ins tabil ity. A second important stabilizing factor, 
discussed in § 15.11 is the expectation that the tower should 
quickly transition to the relativistic regime. Once the outflow 
becomes ultra-relativistic with a very large 7-factor, the rela- 
tivistic time delay may effectively stabilize the tower (Gian- 
nios & Spruit 2006). The reason for this is that MHD instabil- 
ities grow on the local Alfven-crossing time in the fluid frame 
and hence much slower in the laboratory frame. As a result, 
even if instabilities are excited, they do not have enough time 
to develop before the break-out of the flow from the star. Fi- 
nally, the magnetic tower is not in a vacuum, but is being con- 
fined by a high-pressure cocoon. The cocoon may provide 
some stabilization at least for external kink modes, although 
internal kink modes will not be affected. Some theoretical ev- 
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idence supporting the idea of the external pressure stabiliza- 
tion can be derived from Konigl & Choudhuri's (1985) analy- 
sis of a force-free magnetized jet that confined by an external 
pressure. Indeed, they showed that a non-axisymmetric he- 
lical equilibrium state becomes energetically favorable (con- 
serving the total magnetic helicity in the jet) only when the 
pressure drops below a certain critical value. If this hap- 
pens and the external kink does go unstable, then this non- 
axisymmetric equilibrium may perhaps be interpreted as the 
end point of the non-linear development of the instability. 

5.5. Reconection Inside the Tower 

Another important process that may affect the propagation 
of the magnetic tower through the star is magnetic reconnec- 
tion across the cylindrical separatrix current sheet at R = R s . 
This process could, in principle, lead to the break up of a sin- 
gle tower into a train of smaller spheromak-like plasmoids. 
Instead of a further lengthening of the tower, one would then 
effectively get continuous injection of new plasmoids at the 
base of the outflow. The new plasmoids would push the exist- 
ing ones further up the axis, so that the net dynamical effect 
would probably be essentially the same as that of one sin- 
gle continuously growing tower. The physics of this process 
is in fact identical to the cyclic generation of plasmoids via 
reconnection across the separatrix suggested in the context of 
magnetospheres of accreting young stars, resulting in a knotty 
jet (Goodson et al. 1999; see also Uzdensky 2004). At what 
level this tearing instability should saturate, i.e., what should 
be the expected size of the plasmoids, is presently not known. 

We would like to remark, however, that we doubt that a fast 
and efficient reconnection process is possible in the tower, as 
long as it is still deep inside the star. It is indeed true that 
fast reconnection is commonly believed to be the underly- 
ing mechanism for many spectacular and most energetic as- 
trophysical phenomena, from substorms the Earth magneto- 
sphere, to flares in the solar corona (e.g., Sweet 1958), to X- 
ray flares in magnetars (Lyutikov 2003). However, all these 
cases are characterized by relatively low-density environ- 
ments, where classical collisional resistivity is so small that 
other, non-classical terms in the generalized Ohm law (i.e., 
the Hall term or anomalous resistivity) become important. 
This happens when the reconnection layer thickness, com- 
puted from the classical Sweet-Parker theory (Sweet 1958; 
Parker 1957) with Spitzer resistivity, becomes smaller than 
the relevant microscopic plasma length scales, such as the col- 
lisionless ion skin depth and/or the ion gyro-radius. This con- 
dition is easily satisfied both in the solar corona and in the 
Earth magnetosphere. Then, as many numerical simulations 
(e.g., Mandt et al. 1994) have shown, reconnection proceeds 
at a fast rate, in part due to the so-called Petschek mechanism 
(Petschek 1964). 

However, in the dense collapsar environment the situa- 
tion is different. Even if the plasma pressure inside the 
magnetic tower constitutes a tiny but finite fraction of the 
magnetic pressure, the resulting electron-positron pair den- 
sity is so high that the plasma should be considered strongly 
collisional. For example, let us assume, as an illustration, 
that the interior of the tower is filled with an optically-thick 
radiation-pressure dominated pair plasma with the thermal en- 
ergy density equal to 10~ 3 of the magnetic energy density 
and let us make some simple estimates relevant to recon- 
nection physics. For a characteristic magnetic field strength 
of 10 14 G, the corresponding plasma temperature should be 
of order T ~ 3 • 10 9 K ~ 300 keV. The corresponding 



equilibrium pair density [using the non-relativistic formula 
n e = 2{m e kBT /2irh 2 ) 3 / 2 exp(-m e c 2 / ksT) for a rough estimate] 
is on the order of n e <~ 2 • 10 29 cm" 3 . The classical non- 
relativistic Spitzer resistivity due to electron-positron colli- 
sions at the above temperature is r\ <~ 0.1 cm 2 /sec. Now, the 
main dimensionless quantity that characterizes the reconnec- 
tion layer in resistive MHD regime is the Lundquist number 

5 = VaL/tj, where L is the global length-scale of the recon- 
necting system. Taking L = 10 7 cm and using the speed of 
light instead of the Alfven velocity, we get S <~ 3 • 10 18 . Then, 
the classical Sweet-Parker reconnection theory gives us an 
estimate for the thickness of the reconnection layer as small 
as <5sp = LS~ l l 2 <~ 6 • 10~ 3 cm. This thickness is of course 
tiny compared with L, but it is huge compared with the rel- 
evant physical plasma length-scales, such as the electron col- 
lisional mean free path (/ m f p <~ 10" 6 cm), the electron gyro- 
radius (p e ~ 10" 11 cm), and the electron collisionless skin- 
depth (d e = c/ujp e <~ 10~ 9 cm). This means that, although 
they are themselves very small, the classical collisional re- 
sistivity (due to e + - e~ Coulomb collisions) and the Comp- 
ton drag (due to e + - 7 and e~ - 7 collisions) greatly exceed 
the collisionless terms in the generalized Ohm law. There- 
fore, resistive MHD (augmented by the Compton drag which 
may give a comparable, perhaps even somewhat larger than 
Spitzer, contribution to the total resistivity) should provide an 
accurate description of the plasma, even inside a very thin re- 
connection layer. 

During the past two decades it has been convincingly 
demonstrated, via numerical simulations (e.g., Biskamp 1986; 
Uzdensky & Kulsrud 2000), theoretical analysis (Kul- 
srud 2001; Malyshkin et al. 2005), and laboratory plasma 
experiments (Ji et al. 1999), that in collisional resistive-MHD 
systems reconnection proceeds in the very slow Sweet-Parker 
regime and that the Petschek fast reconnection mechanism is 
disabled. If we can extend this result to the highly relativistic 
and optically-thick electron-positron plasma inside the collap- 
sar magnetic tower, we then have to conclude that any large- 
scale reconnection across the separatrix is very slow and in- 
efficient. One has to be aware of several caveats, however. 
First, the role of the photon drag on reconnection dynamics 
is not known, especially its effect on the geometry of recon- 
necting flow. Second, there is a possibility that, even in the 
resistive-MHD regime, reconnection can be greatly enhanced 
in the presence of MHD turbulence, as suggested by Lazarian 

6 Vishniac (1999). And third, since magnetic energy den- 
sity in the tower greatly exceeds both the plasma pressure and 
its rest-mass energy density, reconnection flow is bound to 
be strongly relativistic; then special-relativistic effects, such 
as the Lorentz contraction, may modify the conclusions de- 
rived from the non-relativistic reconnection theory (Blackman 
& Fields 1994; Lyutikov & Uzdensky 2003). With all these 
reservations in mind, we still believe that reconnection pro- 
cesses are strongly inhibited in the tower when it is still deep 
inside the collapsing star. However, as the tower grows and 
eventually breaks out of the star, the plasma cools and the par- 
ticle density in it drops rapidly. Therefore, we should expect 
that the plasma will become collisionless, at least as far as 
reconnection physics is concerned, at late times, perhaps af- 
ter the break-out. This opens up the possibility that any pos- 
sible disruption of the tower through reconnection (and cor- 
responding magnetic energy release) is delayed and happens 
only after the tower clears the stellar surface. 
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5.6. Prospects for Numerical Simulations 

In this paper we have developed a mathematical framework 
describing a sequence of static force-free equilibria which 
represents a magnetic tower inside a star. To verify if our 
picture is realizable in practice and to make further progress 
in understanding the development and evolution of magnetic 
towers in stars, it is desirable to perform time-dependent nu- 
merical simulations. Such simulations will determine the de- 
gree to which the simplifications we have adopted for our an- 
alytic model hold in realistic environments. The key feature 
to be investigated is the structure of the tower and the cocoon 
and how they may be affected by instabilities. 

We envision a sequence a numerical simulations in which 
increasingly realistic physics is incrementally added to the 
simplest elemental simulation. Each of these studies will be 
able to answer a subset of key questions with increasing de- 
gree of realism. First, the following questions can be ad- 
dressed with non-relativistic axisymmetric MHD simulations. 
What basic magnetic field configuration results, in practice, 
when the conditions we describe are set up? Under what con- 
ditions does a magnetic tower form, if at all? How do the 
tower growth and structure change as the tower expands into 
regions of lower density? To what degree is cocoon material 
mixed into the tower, and magnetic field from the tower mixed 
into the cocoon? How strongly is magnetic field concentrated 
toward the axis? How does the Poynting flux depend on ra- 
dius and height? What is the overall field line geometry? Is 
magneto-centrifugal acceleration from the base of the outflow 
inhibited? In simulations where multiple towers are investi- 
gated, how do they interact? To what degree does the cocoon 
help collimate and stabilize the tower? How rapidly do the 
cocoon walls spread laterally? Can the cocoon expansion re- 
sult in the disruption of the star? What differences are there 
between cocoons formed by purely hydrodynamical jets and 
by magnetic towers? 

Next, relativistic MHD simulations will be able to address 
fundamental questions of critical interest to the application of 
magnetic towers to gamma-ray bursts. Of particular interest is 
the beaming angle and angular distribution of energy, Lorentz 
factor and magnetic field after the tower makes the transition 
to relativistic expansion and when it eventually breaks out of 
the stellar surface. Do relativistic effects suppress the growth 
of instabilities in the tower and surrounding plasma? In ad- 
dition, inclusion of relativity is necessary for some regions of 
the central engine, e.g. outside of a magnetar light cylinder or 
the inner parts of black hole accretion disks. 

Finally, three-dimensional simulations, both relativistic and 
non-relativistic, will enable study of non-axisymmetric insta- 
bilities. Intermixing of baryonic material and magnetic field 
between the tower and the cocoon may also be quantitatively 
addressed with three-dimensional simulations. 

Several numerical schemes have recently been developed 
(Koide, Shibata & Kudoh 1999; Gammie, McKinney & Toth 
2003; Del Zanna, Bucciantini & Londrillo 2003; De Vil- 
liers, Hawley & Krolik 2003; Fragile 2005; Komissarov 2005; 
Nishikawa et al. 2005) which should be capable of simulating 
magnetic towers powered by central engines formed after core 
collapse. In fact, general relativistic MHD simulations of col- 
lapsars have already been performed by two groups (Mizuno 
et al. 2004; De Villiers, Staff, & Ouyed 2005). In prin- 
ciple these simulations may contain structures which could 
be termed magnetic towers as described in this paper. How- 
ever, it may be difficult to clearly discern the properties of the 



magnetic structure in complex simulations. For this reason 
we suggest above incremental simulations guided by analytic 
studies as presented here. 

While the magnetic tower itself is force-free, pressure and 
inertia of the surrounding plasma play a key role in its devel- 
opment. This points to the desirability of magnetohydrody- 
namical (MHD) simulations capable of simultaneously solv- 
ing force-free and hydrodynamical regions of the flow. One 
possible choice is a hybrid numerical scheme which solves 
relativistic force-free equations in highly magnetized parts of 
the computational domain, e.g. within the tower, with a mov- 
ing contact discontinuity evolved within a hydrodynamical 
simulation (Blandford 2005, private communication). An- 
other possible strategy is to evolve relativistic MHD equa- 
tions, but switch to relativistic force-free equations whenever 
B 2 becomes larger than pc 2 (e.g., McKinney 2006). 

As an ultimate goal, numerical simulations of the full prob- 
lem, including a detailed description of the central engine with 
relevant microphysical processes and neutrino transport, are 
desirable for a comprehensive understanding of the formation 
and evolution of a magnetic tower in a star. 

6. CONCLUSIONS 

In this paper we have proposed a new mechanism for 
driving a strongly-collimated, baryon-poor, magnetically- 
dominated outflow through a massive star in the collapsar 
model (e.g., MacFadyen & Woosley 1999) for long-duration 
gamma-ray bursts. This magnetic mechanism may work inde- 
pendently and in parallel with the usual neutrino-driven out- 
flow mechanism. It may also be relevant for asymmetric core- 
collapse supernova explosions. 

To model the magnetically-dominated outflow, we invoke 
the concept of the magnetic tower, first introduced by Lynden- 
Bell (1996) in the AGN context. Overall, our model can be 
outlined as follows. The core collapse of a massive rotat- 
ing star may result in two distinct configurations, both plau- 
sible candidates for the GRB central engine. The first one is 
a stellar mass black hole accreting rapidly through an accre- 
tion disk. The second is a millisecond magnetar. A hybrid 
configuration with a disk around a magnetar is also possible. 
In the first case, the core of a rotating massive star collapses 
into a black hole (or forms a proto-neutron star), whereas the 
overlying stellar material that is continuously falling towards 
the center may possess enough angular momentum to form 
an accretion disk. Any pre-existing stellar magnetic field be- 
comes greatly amplified by the combined action of compres- 
sion during the collapse and field-line stretching due to differ- 
ential rotation. In addition, the disk is generally unstable to 
magneto-rotational instability (MRI) and becomes turbulent. 
This results in further amplification of the magnetic field due 
to turbulent dynamo action inside the disk, up to the level that 
may be in excess of 10 16 G (as suggested by Akiyama et al. 
2003 in the core-collapse supernova context). Even though 
this field is predominantly toroidal, the poloidal field compo- 
nent may also be significant; we believe that a figure of 10 15 G 
for the large-scale poloidal magnetic field is not unreasonable. 

As a result of the amplification, magnetic tubes become 
buoyant and escape out of the disk into the overlying low- 
density corona, forming coronal loops. These loops are an- 
chored at both ends at the disk surface but, generally speaking, 
the two footpoints of a loop are at different distances from the 
center and hence rotate at different angular velocities. The re- 
sulting differential rotation further twists the magnetic loops, 
generating more toroidal flux, and the disk magnetosphere 
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tends to inflate. A similar process has also been shown to 
work in the case of a magnetic field linking the disk directly to 
the black hole (Uzdensky 2005). In any case, it is this inflation 
process that is ultimately responsible for the formation of the 
magnetically-dominated outflow in our model. For simplic- 
ity, we consider the evolution of an axisymmetric magnetic 
arcade. Because the density in the corona above the disk is 
low, the initial expansion is force-free and takes place at some 
finite angle with respect to the rotation axis. Subsequently, 
however, the expanding magnetic field weakens to such a de- 
gree that it starts to feel the presence of the ambient stellar 
gas. As was shown by Lynden-Bell (1996), any external gas 
pressure surrounding an expanding twisted magnetic arcade 
eventually becomes dynamically important and channels the 
expansion into the vertical direction; a tall magnetic tower 
forms and continuously grows as a result. In our model, we 
modify Lynden-Bell's model by taking into account that the 
tower expansion is supersonic with respect to the unperturbed 
stellar gas. We thus envision the growing magnetic tower act- 
ing as a piston that drives a strong shock through the star. The 
hot shocked stellar material between the shock and the tower 
forms a high-pressure cocoon, similar to that seen in hydro- 
dynamic simulations by Zhang et al. (2003). This cocoon 
envelopes the tower and provides the collimating pressure for 
it. Thus, in our model, the tower is confined not by the pres- 
sure of the background stellar material, but by its inertia; the 
strong shock and the cocoon act as mediators that convert the 
inertial support into the pressure support that ultimately acts 
on the tower. 

The entire configuration grows vertically with time and 
eventually reaches the star's surface, thus providing a very 
narrow baryon-clean channel in the form of a Poynting-flux 
dominated jet, surrounded by a less-collimated hot gas out- 
flow. One of the great advantages of the magnetic tower 
model in the GRB context is that the magnetic field lines re- 
main closed (e.g., with both ends tied to the disk) during the 
expansion process and, as a result, the volume occupied by 
the tower remains insulated from the surrounding stellar gas 
and hence relatively baryon-free. As the tower grows, it just 
pushes the stellar gas out of its way into the cocoon. 

In addition, we believe that the physical scenario developed 
in this paper may also be applicable to the case where the 
central engine is a rapidly-rotating (i.e., millisecond) magne- 
tar. The rotational energy of the magnetar is extracted mag- 
netically as a Poynting flux that inflates a baryon-free cav- 
ity inside the progenitor star. This mechanism is similar to 
that proposed by Ostriker & Gunn (1971) for powering type II 
supernovae. However, whereas their model was spherically- 
symmetric, we suggest that a relativistic generalization of the 
magnetic tower model can naturally result in splitting of the 
magnetically-dominated outflow into a pair of oppositely di- 
rected narrow jets. 

In our paper we first give a simple description of Lynden- 
Bell's mag netic tower model in its original accretion disk con- 
text (§ 12. Q . Then we present our model of a m agne tic tower 
propagating and driving a shock through a star (§ I2.2I>. We de- 
rive some basic scalings for the tower pa rame ters (§ 12. 3> and 
make some simple estimates for them f§ 12.41 . Then, in §|3]- 
|4j we illustrate our model by two specific analytical solutions 
for the magnetic field s truc ture in the main body of the tower; 
one of the solutions (§ I4.2> is characterized by a very narrow 
axial jet with a non-zero poloidal line current present in the 
core of the tower. In §13 we briefly discuss several interest- 
ing physical issues that, although important, lie beyond the 



scope of the present paper. Namely, in § 15.11 we argue that, 
even though it may start off non-relativistically, the expansion 
of the tower will quickly have to transition to the relativistic 
regime. In § 15.31 we discuss the issue of bary on contamina- 
tion of the tower by disk winds. Next, in § 15 .41 we touch upon 
the question of stability of the magnetic tower configuration 
and address the ro le th at ideal MHD instabilities may play in 
our scenario. In § !5.5l we discuss the possibility of the tower 
disruption and its break-up into a series of smaller plasmoids 
via internal reconnection. We argue, however, that, while the 
tower is still deep inside the star, the plasma inside of it is 
likely to be strongly collisional, so that classical resistive- 
MHD should apply inside the reconnection region. If this is 
the case, then, we believe, any reconnection process should 
be very slow, and so a global disruption through reconnection 
is not likely during this stage. During subsequent expansion, 
however, the density drops and the reconnection process may 
enter a faster collisionless regime, perhaps resulting in rapid 
magnetic energy dissipation. 

Finally, in accordance with our discussion in § [5] we can 
formulate the following three areas that, we believe, require 
immediate attention and should and can be addressed in the 
near future. They are: (1) a special-relativistic generaliza- 
tion of the magnetic tower model; (2) 3D relativistic MHD 
numerical simulations of a magnetic tower inside a star; (3) 
the conditions for and the non-linear development of, various 
MHD instabilities. 
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